Energy expressions for n = 3 and 4 systems in a single-j shell 
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For systems with three and four fermions within a single-j shell, analytical expressions for the 
state energies are presented from a decomposition of the angular momentum. In some important 
cases the expressions acquire a very simple form. The expression may help us in understanding the 
structure of isomeric states. The decomposition also makes it possible to construct the algebraic 
condition for conservation of seniority. 
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The developments in experimental techniques have 
made it possible to study the structure oi N Z nuclei 
far from the stability line. For N k, Z nuclei just above 
''"Ca and immediately below ^°°Sn, the structure of most 
low-lying states is dominated by the coupling of few va- 
lence particles (or holes) moving in a single-j shell-model 
orbit. In a single-j shell, a state with angular momentum 
/ may be constructed for which the expectation value of 
a two-body interaction can be written as. 



Ei = C'jVj, 



(1) 



where Vj = (j^; J|V^|j^; J) are two-body matrix elements 
and Cj are corresponding expansion coefficients (The 
constant contribution from the single particle energy has 
been neglected for simplicity). If isopsin symmetry is 
conserved in the two-body interaction V, then one has a 
total number of 2j -I- 1 matrix elements being coupled to 
spin values J ^ to 2j [U, Q ■ The expansion coefficient 
C^gives the number of pairs with angular momentum 
J [3| . Some trivial examples are the two-nucleon system, 
for which Ej = Vj, and the case of a closed single-j shell, 
for which Ej = Y,j{2J + l)Vj [i!, 0. 

The total number of pairs with all spins J is given 
by Si, 



and 



J,odd 



\ [n 
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(3) 



where n is the total number of nucleons and T is the total 
isospin quantum number of the system. 

In this work systems with n = 3 and 4 nucleons are 
analyzed and expressions for the number of pairs with 
any total angular momentum and isospin are derived in 
terms of angular momentum coupling coefficients. Sim- 
plified relations are also derived when possible. For ear- 
lier works based on coefficients of fractional parentage 



calculations, see Refs. 0-13. In Ref. [31], the number of 
J = pairs in the 0+ states of the n = 4 nucleus '^"'Ti is 
obtained. We will start from the so-called multistep shell 
model Q. It employs correlated two-particle bases as 
building blocks and was first proposed as a truncation of 
the exact shell-model approach [y] (The generalization of 
this model to many-body systems is given in Refs. 
When restricted to a single-j shell, the correlated basis 
with spin J is closely related to the corresponding ma- 
trix element Vj, which makes the multistep shell model 
a very convenient tool to study the energy expressions of 
nuclear systems. 

For a n = 4 system with two protons and two neutrons, 
the system can be decomposed into proton and neutron 
blocks. The wave function of a given state with total 
angular momentum / can be written as Q, 



(4) 
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where Xj(JpJn) is the amplitude of the four-body wave 
function and Jp and J„ are even numbers denoting the 
angular momenta of the proton and neutron pairs, re- 
spectively. 

The four nucleons can couple to spin / = to 2(2j — 1) 
and isospin T ~ 0, 1 and 2. The T ~ 2 state is the double 
analog of the system of four identical nucleons. 

The dimension of a given state / will be denoted by 
D/(j), while for the number of states with a given isospin 
T the notation Dj{j) will be used. Recently, there have 
been efforts to obtain algebraic formulas for the state 
dimension in a single-j shell p^| - [l3 |. For / = and 
n = 4 these quantities are relatively easy to evaluate. 
The dimension equals to the number of T = 1 states in 
the two-body system, i.e., Do{j) = (2j + l)/2. For spin 
zero states with T = 2 one has D^{j) = [(2j + 3)/6] = 
[(j -I- 2)/3] {[n] denotes the largest integer not exceeding 
n) [13, m. For the T = states it is D^{j) = Do{j) ~ 
Dl{j) = (2j-hl)/2- [(2j-h3)/6]. That is, for j = 3/2 to 
7/2 one gets I?g(j) = (2j - l)/2 and Dl{j) = 1. For j = 
9/2 to 13/2, we have D[J(j) = (2j - 3)/2 and Dl{j) = 2. 
For / = 0, there is no T = 1 state. 

In the bases of Eq. (|4]), the matrix elements of the 
Hamiltonian can be expressed as linear combinations of 



two-body interactions Vj as, 

{jUJp)jl{Jn);i\v\jl{4)fM;i) 

= {Vj, + VjJSj^j,Sj^j^ + ^ A/j(Jp J„; j; j;)K7, (5) 

J 

where the spin J can take both even and odd values (J ~ 
to 2j). The matrix M is given as 
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where j = v2J+T and A and j are half integers. Both 
the Hamiltonian matrix and M are symmetric since we 
have Mj{JpJn; J'^J'^-^) = Mj{J^j;^; JpJn). 

Direct diagonalization of the Hamiltonian matrix 
(Eq. ([5])) would give cigen energies as nonlinear combi- 
nations of two-body interactions Vj. But still the eigen 
energy can be re-expressed in the linear expansion form of 
Eq. (1). Generally, the expansion coefficients for a state 
with total angular momentum / of the n = 4 system can 
be written as. 



■Jr,Jn:J'J' 



where the first and second terms in the right-hand 
side give contributions from identical nucleon pairs and 
proton-neutron pairs, respectively. The wave function 
amplitudes X as well as the coefficients C can be cal- 
culated numerically by diagonalizing the Hamiltonian 
matrix. The thus calculated coefficients are usually ir- 
rational numbers and depend on the interactions em- 
ployed 

With the expression given above, calculations have 
been done to explore the energy expressions in several 
single-j shells. As examples, in the upper panel of Fig. [1] 
coefficients Cj are plotted as a function of spin / for yrast 
states of ^^Cd. Calculations are done in the Ogo^ shell 
with two-body interactions extracted from Ref. [i5| . 

Eq. ([7]) is valid for all spins and isospins. But in some 
cases, the expression can be greatly simplified based on 
the symmetric properties of the 6j symbols. It is easily 
seen that for J = we have rank(M) = 1 and, 

Af°(JpJ„; J;j:J = aj{Jp)aj{j;)Sj^jJj,j^, (8) 




FIG. 1: (Color online) Upper: expansion coefRcients Cj as a 
function of spin I for yrast states of ^®Cd in the OgQ/2 shell; 
Lower: same as the upper panel but for those of '^'^Sc (and 
■*^Ti) in the O/7/2 shell. 



For the number of J = pairs in system with any angular 
momentum / one gets. 



' 7' /*'2 J Jp Jn I 



MQ{JpJn', JpJn) — '^JpJnJpJn/j 



3 3] 
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Jp J„ 1 

j j 3 



(10) 



For J = and / = 0, the expression can be even simpler. 
It is ao(Jp) 2Jp/j- and 



Co" = 2|Xo(00)|^ 



JpJpXo{JpJp)Xo{JpJp). (11) 



Another interesting case is the state with the maximum 
spin of / = 4j — 2 and dimension D = 1. For this unique 
state we have Jp = J„ = 2j — 1 and |X/(JpJ„)| = 1. 
Inserting these into Eq. ((T]), it can be easily seen that we 
have Cj = and 3 for the number of pairs with J < 2j— 2 
and J = 2j — 1, respectively. For that of J = 2j — 2 we 
have. 



where 
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C 



I=4j-2 
J=2j-2 



12(4.7- 1)2(4.7 -3)(2j-l) 

2j-l 2j-l Aj-2 
3j - 2 J j 

3j - 2 J 23 - 1 



(12) 
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The number of spin 2j pairs is determined as C'j^jj ^ ~ 

3 — C'j~2j-2- These numbers of pairs are independent of 
interactions used. 

The 12+ state in ^^Fe (four-hole system) has been ob- 
served to be a spin trap [l6| , whose excitation energy is 
lower than that of 10+. This scheme can be reproduced 
nicely in large-scale fp shell model calculations with re- 
alistic nucleon-nucleon interactions [l3]- In O/7/2 shell, 
the energy of this unique state is given as, 



6 



i2( Fc) = ^1/5 + 31/6 



33 
13 



V7, 



(13) 



where V denotes hole-hole interaction matrix elements 
which can be extracted from the spectrum of ^'^Co 
We have a total number of three 10+ states (among which 
one with T = 1). The energy of the first one is given as, 



^10+ (''Fc) 



O.3IOF3 + 1.429^4 
-f 0.497t4i + 1.571V6 + 2.193F7. (14) 

which is calculated to be lower than that of the 12+ state. 
As a comparison, with two-b ody interactions extracted 
from the spectrum of '^^Sc 0, [HI, the energy of the first 
10+ state in ^'*Ti is calculated to be, 

£^io+C''*Ti) = 0.245^3 + 1-473^4 

-fG.570V5 + 1.526^6 + 2.1861/7- (15) 

The expansion coefficients of above two expressions sim- 
ilar to each other, which are obtained by diagonalizing 
the corresponding Hamiltonian matrices. 

A similar spin trap, with I'^ = 16+, has been predicted 
in the four- hole system of ^^Cd below the doubly magic 
^™Sn [l^. This prediction is supported by our Off3/2-shell 
calculations with hole-hole interactions from Ref. [15| as 
well as large-space calculations in the /pg-shcU. In j = 
9/2 shell, the energy of this state is. 



43. 



EierCd) = —V7 + 3Vs + —Vs. 



(16) 



We have a total number of three 14+ states. The energy 
of the first 14+ state is calculated to be, 

E^^+C^Cd) = 0.3971/5 + 1-428V6 

H-0.4931/7 -t- 1.5721/8 -t- 2.200%, (17) 

For the matrix elements we have I/5 w Vy and Vg ~ 
Vs [H, [i^]. Thus the position of the 16+ state relative 
to the first 14+ state is sensitive to the strength of inter- 
action matrix element Vg which is much more attractive 
than above mentioned elements. This observation has 
also been found in Ref. [l9j with numerical calculations. 

For n = 3 and jT^j = 1/2 systems in a single-j shell, 
similar to Eq. Q , the wave function of a state with total 
angular momentum / can be written as, 



Jc 



XiiJ^)\f{J^)j;I), 



(18) 



where Xi{Ja) is the amplitude and Ja and / are integers 
and half integers, respectively. In above equation, we 
have assumed that the system was decomposed into a 
nucleon pair and a odd nucleon. The Hamiltonian matrix 
can be constructed in a similar way to Eq. ([5]). 

In the system with n = 3, one has a total number of 
three pairs with different spins. For the number of pairs 
with a given J one gets. 



c5 



\XjiJ^)\Hj^j 

+ Xj{J^)N'j{J^fjXj{j'J. 



The symmetric matrix N is given as. 



j I J 

J J JL 



(19) 



(20) 



As an illustration, in the lowe panel of Fig. [T] we gave 
calculations for coefficients Cj along the yrast states of 
^•^Sc and ^'^Ti. Calculations are done in the O/7/2 shell 
with two-body interactions from Ref. 

For I j, there is no J = pair. For I = j one 
has NQ[JaJ'a) = 2JaJa/j^ and the expression for the 
number of J = pairs can be simplified as, 



^0 



l^.-(0)l' + ^ E jJLX,iJc.)X,iJ'J. (21) 
^ J.J' 



For the unique state with the maximum spin of / = 3j — 1, 
it can be easily seen from Eq. (|T9)) that we have Cj — 
for J < 2j — 1 and 



C 



1=33-1 

J=2]-l 



r<J=3o-i 



(22) 



The states with total spins / = 1/2 and / = 3j — 2 are 
also unique. For the I ~ 3j — 2 state we have Cj ~ 
and 3/2 for J < 2j — 2 and J = 2j — 1, respectively. For 
the numbers of J ~ 2j — 2 and 2j pairs, they are related 
to the 6j symbol as. 



= 2(4j-l)(2J+l) 



J 3.7-2 J 
j 3 2j - 1 



(23) 

As an example, the energy of the / — 17/2 state in j = 
7/2 shell can be expanded as, 

19 3 in 

E,r/2{J = 7/2) = -V, + -V, + -Vr. (24) 

For the / = 1/2 state, the energy can be expressed as, 

3 



El/2{]) 



yj-1/2 



j + 1/2- 



(25) 



The energy expressions for systems with three and four 
identical nucleons in a single-j shell have been explored 
in textbooks as well as in recent publications [2TI - 
[26j . For four identical nucleons in a single-j shell, the 
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basis state can still be written as \j'^{Ja)j'^{Jp)', I) where 
it is assumed that the system was decomposed into two 
blocks. Since the Pauli principle is not acting, the overlap 
between such states is [y], 

A^ij^jp-^f^ff,) = {f{j^)fij,y,i\fi.rjfi.r^);i) 

J J J Oi. 

-AjJpJ'J'p{ j j ),{2&) 



T' T' T 



where Ja and arc even numbers. For / 
tains, 



one ob- 



A^o(^o J.; JM - 2Jj„ + 4J„ J4 (27) 

where A^(00;00) = (4j - 2)/{2j + 1). In orbits with 
j < 7/2, the dimension of spin zero states is D^{i) ~ 1. 
We found that the energy of this state can be written as. 



Eo{j)= J2 AliJJ;JJ)Vj. 



(28) 



J.cven 



Similarly, the basis states of systems with three iden- 
tical nucleon can be given as the coup ling of a nucleon 
pair and a single nucleon as \j^{Ja)j] I) [231 • These bases 
are not orthonormal to each other. The overlap between 
different states is given as. 



B'j{J.,J'J = {f{Ja)j;i\f{J'Jr:i) 



= ^.JcJL + 2 Jq 



(29) 



It is easily seen that for I ^ j one has Bj{JaJ'g,) = 
1 / 2Aq{J aJ aJ'aJ'a) ■ Similar to systems with four identical 
nucleons, for J < 7/2 we have. 



Ej{:i)=Y,B]{J^J^)Vj, 



(30) 



The dimension oi J = j states in three-identical- 
nucleon system equals to the dimension of J = states 
in four-identical-nucleon system, which is given as = 
[{2j + 3)/6]. It equals to the rank of overlap matrices 
A and B. The dimensions of matrices A and B are the 
same, i.e., n = (2j + l)/2. 

Seniority remains a good quantum number for systems 
with j < 7/2 p. For 9/2 < j < 13/2, there arc two 
I = j and / — states in thrce-identical-nucleon and 
four-identical-nucleon systems, respectively. In this case 
for interactions to conserve seniority, it has to satisfy one 
(or [{2j - 3)/6]) linearly relationship P, [H, [H. In a 
n = 3 system, the conservation of seniority implies that 



{v = 1; j\V\i^ = 3; j)=0 [l|. Based on the overlap matrix, 
we found that the conservation condition can be written 



Y,Bj{OJ)iBj{OJ) 



Bj(OO) 



BU2J)]Vj = 0. (31) 



V 5i(20) ' J 

For interaction Vq, the corresponding coefficient is zero. 
The same conservation condition can be got by employ- 
ing the overlap matrix A of the four-body system. This 
relation is sufficient to ensure that seniority is conserved 
for all states in shells with 9/2 < j < 13/2. If seniority 
is conserved in the interaction, we found that the wave 
function of unique v = 1, I = j state for three-identical- 
nucleon system in any single-j shell satisfies the relation. 



BjiOJa) 



(32) 



It implies that the wave function can be written as, 

l^Ai^ = 1)) = ^^=^i?j(0J„)|j2(J„)j;j). (33) 
Sj(OO) j„ 

The wave function of the = 0, / = state in four- 
identical-nucleon system can be constructed in the same 
way. 

In summary, in this work the number of nucleon pairs 
for three- and four-particle systems in a single-j shell 
has been evaluated with the decomposition of the angu- 
lar momentum in the framework of the multistep shell 
model. We have derived analytic expressions that are 
valid for all spins and isospins. Possible simplified re- 
lations are also explored with the symmetric properties 
of the angular momentum coupling coefficients. With 
the expressions thus derived we explored the number of 
pairs in isomeric states as well as other unique (dimen- 
sion one) states. The decomposition can also help us 
in understanding the conservation of seniority in single-j 
shells. 

Within the framework of the multistep shell model , 
where the building blocks are different nucleon pairs, this 
analysis of pair numbers can be extended to systems with 
many particles. These may help us to understand the 
property oi N Z nuclei in single-j shells as well as 
the structure of isomeric states in many-particle nuclear 
systems, e.g., the spin-trap isomers in ^^Fe and ^^Pd. 
Further work in this direction is underway. 
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The author is supported by a grant from the Swedish Re- 
search Council (VR). 
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